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[. INTRODUCTION 



A. DEVELOPING A MATHEMATICAL MODEL FOR OPTIMIZATION 

Design optimization is a powerful tool available to the design engineer. While 
it is not a new concept, only recently has numerical design optimization become a 
technique which is routinely applied by practicing engineers to the myriad of design 
tasks. This trend can be attributed to recent technological advances in high-speed 
digital computers, which are necessary to solve the multitude of equations modeling 
the optimization problem. 

One area particularly suited for design optimization applications is structural 
design. The design of a structure requires a great deal of forethought by the design 
engineer. Several restrictions (constraints) on the design, such as environmental 
effects, geometry, and material selection, must first be determined. Then, these 
restrictions are applied to the design process to produce a specific output (objective), 
such as maximum structural response or minimum cost. The goal of the design 
engineer is to design the best possible structure, in terms of a specific desired output, 
while adhering to several expected constraints on the final design. That is, the final 
structure must be the optimum design for anticipated conditions. 

In order to achieve this optimum design, the design optimization process 
requires that the design task be defined as a mathematical model. This model 
consists of an objective function, constraint functions, and side constraint functions. 
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These three functions are mathematical equations expressed in terms of design 
variables and state variables. State variables are fixed quantities in the mathematical 
equations. Design variables are those quantities which are allowed to vary during 
the optimization process. The objective function, which contains every design 
variable, is a single equation representing the quantity to be optimized (such as 
weight, displacement, cost, etc.). Constraint functions are one or more equations 
which restrict the design variables; one or more design variables are contained in 
each constraint function. Side constraint functions are equations which limit the 
upper and lower bound ranges of each design variable. Simply stated, the model 
describes, mathematically, what is to be optimized (objective) and the limitations 
(constraints and side constraints) on design variables. 

Having converted the design task into a mathematical model, the design 
engineer can program the model. Using his own optimization code, or, more likely, 
available software, the design engineer may now execute the process of optimizing 
the objective function and the associated design variables. 

B. MODELING A CYLINDRICAL SHELL FOR MINIMUM WEIGHT DESIGN 

The design optimization task is to design a minimum weight thin cylindrical 
shell subjected to external hydrostatic pressure. The shell can be either monocoque 
or ring-stiffened (by internal rectangular cross-section rings placed at equidistant 
intervals within the shell). A single isotropic shell and ring material will be used in 
either case. A more complex orthotropic material or a hybrid construction material 
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could have been used (and supported by the DAPS3 analysis mode). However, due 
to time limitations, only an isotropic material will be investigated. This simpler 
material should still credibly serve the purpose of demonstrating the application of 
numerical optimization techniques to a minimum weight shell design. The optimized 
design is also constrained by buckling, strength, and geometry considerations. This 
design task is applicable to the design of submerged cylindrical structures under a 
static pressure load. Examples might include submarine pressure hulls and torpedo 
cases. 

In order to model the above design task, a shell analysis program, DAPS3 
(Design and Analysis of Plastic Shells version 3), was used to evaluate the buckling 
and strength constraint functions of a thin shell under external hydrostatic pressure. 
Development of these constraint functions is discussed in Chapter II (Shell Buckling 
and Strength). The optimization problem was evaluated by the general purpose 
numerical optimization program ADS (Automated Design Synthesis version 1.10). 
Development of the mathematical model for optimization is discussed in Chapter III 
(Optimization Problem). Both programs, which are written in the FORTRAN 
language, are incorporated into a single FORTRAN code named THESIS. 

The THESIS code produces a minimum weight shell design, and it indicates the 
shell dimension design variables (thickness, ring height, and ring width) associated 
with that optimum design. Therefore, the design task is modeled and optimized 
within the THESIS code. 
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C. THESIS OBJECTIVES 



DAPS3 uses an iterative technique to optimize the shell design variables. The 
variables are manipulated (in a "DO" loop sense) until a design is produced which 
does not buckle under the input hydrostatic pressure. When all possible buckling 
loads are greater than the hydrostatic pressure, DAPS3 terminates the iterations. 
The DAPS3 generated optimum design is produced from the last design iteration. 
It was felt that incorporation of the more sophisticated numerical optimization 
technique should improve upon the DAPS3 iterative design technique. 

In pursuit of the above hypothesis, this thesis will investigate optimum designs 
generated by the THESIS code over a wide range of geometry configurations (L/OD 
ratios) and loads (external hydrostatic pressures). These designs will then be 
compared to the optimum designs generated by DAPS3 (when used in the "design" 
mode). Specifically, the following design elements will be compared: 1) objective 
(shell weight), 2) design variables (thickness, ring height, and ring width), and 3) 
constraints (buckling, strength, and geometry). Any design improvements made by 
use of the numerical optimization technique will be identified. Results are presented 
in Chapter IV. Conclusions as to the advantages and disadvantages of incorporating 
a numerical optimization technique are presented in Chapter V. 



II. SHELL BUCKLING AND STRENGTH 



The purpose of this chapter is to explain the formulation of the buckling and 
strength (or stress) terms used by the THESIS code. A detailed derivation of the 
buckling equations is provided by Renzi [Refs. 1,2]. 

A. BUCKLING EQUATIONS 

Structural failure by buckling is associated with an "unstable" design. Material, 
geometry, and load factors contribute to this instability. The THESIS code designs 
the geometry of a "stable" cylindrical shell based on material and load factors which 
are input by the user. Within the THESIS code, DAPS3 is used as a subroutine 
analysis program to calculate the buckling loads on either a monocoque or ring- 
stiffened shell. These buckling loads are then returned to the numerical optimization 
program (ADS), which is incorporated in the THESIS code, to be considered in 
buckling constraints on the optimum shell design. Of course, other constraints 
(strength and geometry) are also considered before the THESIS code produces the 
optimum shell design. 

The three types of shell buckling loads considered are: 

1. Axisymmetric Collapse Pressure (PC) 

2. Elastic Monocoque Shell Instability Pressure or Elastic Interbay Instability 
Pressure (PIB) 
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3. Elastic General Instability Pressure (PGEN) 



Note that the above buckling terms are expressed in terms of pressure rather than 
in terms of force. The specific unit of pressure used in the THESIS code is (lb/in 2 ). 
The use of buckling pressures, rather than buckling forces, allows easy comparison 
of the buckling loads to the applied load, which is the external hydrostatic pressure 
(PSTR), also in (lb f /in 2 ). The comparison of buckling load to applied load will 
quickly reveal whether the shell will buckle under the specified applied load. From 
herein, any use of the term "load" will mean pressure in units of (lb f /in 2 ). 

Since this thesis examines both monocoque and ring-stiffened shells, the type 
of shell also determines which particular buckling loads are present. For the 
monocoque shell and ring-stiffened shells, two buckling loads are present: 
axisymmetric collapse pressure (PC) and elastic shell monocoque shell instability 
pressure (PIB). When the shell is ring-stiffened, a third buckling load, elastic general 
instability pressure (PGEN), is also present. 

Note that the buckling term PIB is associated with two different names. If the 
shell is monocoque, then PIB is called the elastic shell monocoque instability 
pressure. Otherwise, the same term is called the elastic interbay instability pressure 
if the shell is ring-stiffened. This distinction is discussed in the subsection concerning 
development of the PIB term. 
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1. Axi symmetric Collapse Pressure (PC) 

This buckling load is present in both monocoque and ring-stiffened shells. 
It describes local axi symmetric buckling of the shell in between two ring-stiffeners 
(or in between the end frames for a monocoque shell). An example of axi symmetric 
buckling between two rings is illustrated in Figure 1. 




Figure 1. Local Axisymmetric Buckling of a Ring-Stiffened Shell Between 

Adjacent Rings. 

To find an explicit expression for the axisymmetric collapse pressure (PC), 
first consider the differential element of a thin cylindrical shell shown in Figure 2. 
After establishing equilibrium of forces and moments on the element, it is found that 
the governing differential equation of the isotropic shell is: 
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where: 

T = Thickness of Shell 

R Q 1 1 

a = — ; R c = —(Outside Diameter), R - R c (7) 

R 2 2 



v = Poisson Ratio 



E = Elastic Modulus 
w = Radial Displacement 




Figure 2. Differential Element of an Isotropic Thin Shell 

Subjected to External Hydrostatic Pressure (PSTR). 



The coefficient of the second order term, in this fourth order ordinary 
differential equation, renders the solution w(x) to be nonlinear with respect to the 



8 



